Evidence for a direct band gap in the topological insulator Bi2Se3 from theory and 

experiment 
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Using angle-resolved photoelectron spectroscopy and ab-initio GW calculations, we unambigu- 
ously show that the widely investigated three-dimensional topological insulator Bi2Se3 has a direct 
band gap at the V point. Experimentally, this is shown by a three-dimensional band mapping in 
large fractions of the Brillouin zone. Theoretically, we demonstrate that the valence band maxi- 
mum is located at the Brillouin center only if many-body effects are included in the calculation. 
Otherwise, it is found in a high-symmetry mirror plane away from the zone center. 

PACS numbers: 71.15.m, 71.20.b, 71.70.Ej, 79.60.-i 



Bismuth selenide is a prototypical topological insula- 
tor (TI) that, after its identification as such a material 
[l|,[2[, has been very widely studied. Its surface electronic 
structure consists of a single Dirac cone around the sur- 
face Brillouin zone (SBZ) centre f , with the Dirac point 
(DP) placed closely above the bulk valence band states. 
In order to exploit the multitude of interesting phenom- 
ena associated with the topological surface states Q , 
it is necessary to access the topological transport regime 
with the chemical potential near the DP, but simultane- 
ously in the absolute bulk band gap. Due to the close 
proximity of the DP and the bulk valence states at f , 
this is only possible if there are no other valence states in 
Bi2Se3 with energies close to or higher than the DP. As 
the bulk conduction band minimum (CBM) is undisput- 
edly placed at F (see, e.g., Refs. ^ and'o'), the question if 
the absolute valence band maximum is below the DP at 
F or at other locations in the Brillouin zone (BZ) is iden- 
tical to the question about the nature of the fundamental 
band gap, direct or indirect. 

The nature of the bulk band gap is thus of crucial 
importance for the possibility to exploit the topologi- 
cal surface states in transport but the position of the 
VBM in band structure calculations remains disputed. 
In a linearized muffin-tin orbital method (LMTO) calcu- 
lation within the local density approximation (LDA), the 
VBM was found at the F point, implying that Bi2Se3 is 
a direct-gap semiconductor 0. Contrarily, by employ- 
ing the full-potential linearized augmented-plane-wave 



method (FLAPW) within the generalized gradient ap- 
proximation (GGA), the authors of Ref. 5 have found 
the VBM to be located on the Z - F line of the BZ, 
which is lying in the mirror plane. Similar results have 
been obtained in Ref. [sl with the plane- wave pseudopo- 
tential method (PWP) within the LDA. Various density 
functional theory (DFT) calculations of the surface band 
structure of Bi2Se3 [l|, IH, IH, |lo| also indicate that the 
VBM of bulk bismuth selenide is not located at the BZ 
center. In these calculations the projected bulk states 
show the VBM to lie on the F - M direction of the SBZ. 
A shift of the VBM towards the F point caused by tak- 
ing into account many-body corrections within the GW 
approximation was demonstrated in Ref. 8. So far, how- 
ever, the valence band dispersion has only been calcu- 
lated along high-symmetry lines of the BZ, precluding a 
firm conclusion about the position of the VBM in the 
three-dimensional BZ. Thus, the nature of the band gap 
has not been determined unequivocally so far. 

Experimentally, it appears that the assumption that 
bismuth selenide is a material with a direct band gap has 
first been made in Ref. [ll|, where galvanomagnetic prop- 
erties were investigated. By inspecting angle-resolved 
photoemission spectra in the T — Z — F direction, the 
authors of Ref. |5| have concluded that the VBM is lo- 
cated at the F point. More recent experimental evidence 
appears to support this but most studies have focused 
on the topological state in the immediate vicinity of F 
(see e.g. Ref. lil). Even if a systematic exploration of 
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FIG. 1: (Color online) Photoemission intensity and theoretical results at different constant parallel momenta /cf^ along 
directions parallel to V — Z (a) and normal momenta kz along directions parallel to Z — U (b), which are shown by vertical 
and horizontal dashed lines in Fig. O respectively. Dashed horizontal lines correspond to the experimentally observed energies 
of the VBM (0.505 db 0.030 eV) and the CBM (0.170 ± 0.025 eV). The theoretical curves (shifted to have the VBM at the 
same energy as in the ARPES experiment) present the lowest conduction band and the uppermost valence band obtained for 
the experimental [13] (unrelaxed) and relaxed atomic positions without (LDA) and with (LDA+GVK or GGA+GVK) the GW 
corrections to the DFT (LDA or GGA) band structure. The two parallel and nearly vertical lines of higher photoemission 
intensity seen at kz values of roughly 2.8 and 3.0 A in Fig. (a) correspond to Bi 5(i core levels which have been excited by 
second-order light from the monochromator. 



the experimental dispersion along the F — M is lacking, 
the existing evidence thus points towards a VBM at the 
F point, in contrast to most first-principle calculations 
that place the VBM quite far away from the F-point. 

In this Letter, we report on theoretical and experimen- 
tal evidence for a direct band gap in bismuth selenide. 
We systematically explore the valence band structure in 
large fractions of the BZ and show theoretically that the 
VBM is found at F, but only if GW corrections are taken 
into account. The calculated dispersion for the upper VB 
is found to be in good agreement with ARPES results. 
Our findings are important for the scattering and trans- 
port properties of the topological surface states and may 
call for a new interpretation of previous results in. e.g. 
quasiparticle interference from these states, which involve 
the bulk band-structure features derived within the DFT 
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Experiments were performed on single crystals of 



Bi2Se3 [1^. ARPES measurements were carried out on 



the SGM-3 beamline of the ASTRID synchrotron radi- 
ation facility [l6|. In order to probe the Bi2Se3 bulk 
band structure, the photon energy was varied from 14 to 
32 eV in 0.1 eV steps. The combined energy resolution 
was < 22 meV and the angular resolution was ~0.15°. 
Samples were cleaved in-situ at room temperature, then 
cooled to ~70 K for measurements. The kz values for 
the three-dimensional representation of the photon en- 
ergy scan were calculated using free-electron final states. 



i.e., kz = \j2me/h'^{Vo + ^^j^^^^ cos^ (6>))^/^, where is 
the emission angle and Vq is the inner potential. From 
the normal-emission ARPES data it is easy to locate a F 
point. The value of Vb = 11.8 eV was determined by iter- 
atively changing Vq such that the kz value corresponding 
to this F point agrees with the size of the Bi2Se3 BZ. It 
should be noted that = 11-8 eV is in good agreement 
with values reported previously 13, [H, [l3, [isl and is the 
only value between 1 and 26 eV which results in a correct 
kz value for the F point. 

Experimentally, we go beyond the standard approach 
of ARPES band structure determinations by probing a 
large fraction of /c-space on a dense grid of emission an- 
gles and photon energies, even though only a fraction of 
the data in the most relevant f — M direction of the SBZ 
is shown here. Combined with the assumption of free 
electron final states, this permits us to plot the photoe- 
mission intensity as a function of kz for a given /cy = /cf ^ 
value parallel to the surface along F — M [Fig.[T]^a)] or at 
constant kz as a function of /cf^ [Fig.fljb)]. This uncon- 
ventional way of plotting ARPES results is only possible 
based on a large data set, but it is excellently suited for 
a comparison to band structure calculations. 

Our ab initio calculations were performed by employ- 
ing the FLAPW method as implemented in the FLEUR 
code[iiwithin both the LDA of Ref. H and the GGA 
of Ref. '21I for the exchange-correlation (XC) functional. 
The mentioned two approximations (LDA and GGA) 
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FIG. 2: (Color online) Contour plots of the lowest conduc- 
tion band (upper row) and the uppermost valence band (lower 
row) in the mirror plane (see Fig. SI of Supplemental Material 
[23I]). The presented results are obtained for the experimen- 
tal atomic positions without (LDA) and with (LDA+CVK) 
the GW corrections to the LDA band structure. Vertical 
and horizontal dashed lines drawn on the valence-band con- 
tour plot in the LDA-\-GW case correspond to the cuts of the 
three-dimensional photon-energy scan, which are presented in 

Fig.m 



have been used in order to reveal the effect of different 
reference one-particle band structures on the GW results. 
The GW approximation (with the inclusion of spin-orbit 



interaction as implemented in the SPEX code [13|, |22|) 
was applied in the one-shot framework, where the Kohn- 
Sham eigenfunctions are taken as approximate quasipar- 
ticle wave functions. As an additional factor that can 
affect our ab initio results, we consider different atomic 
positions for Bi and Se atoms in a rhombohedral crys- 
tal structure with experimental lattice parameters taken 
from Ref. 
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One set of positions corresponds to the 
atomic positions reported in Ref. fl2l (labeled as "unre- 
laxed"). Another set (labeled as "relaxed") was obtained 
during a relaxation procedure optimizing the atomic po- 
sitions at fixed volume. For more computational details, 
we refer the reader to Sec. SI of the Supplemental Mate- 
rial 

It should be noted that our GW study goes beyond 
the one of Ref. H in several aspects: Instead of the PWP, 
we employ the FLAPW method, which treats core, va- 
lence, and conduction electrons on an equal footing. Fur- 
thermore, we do not resort to a plasmon-pole model for 
the dielectric matrix but use the random-phase approxi- 
mation (RPA) explicitly without approximations for the 
frequency dependence. In Ref. the GW calculation was 



performed without spin-orbit coupling; spin-orbit inter- 
action was included by second variation using the GW 
corrected eigenvalues only after the quasiparticle spec- 
trum had been obtained. We use the full four-component 
spinor wave functions, as obtained from a fully relativis- 
tic DFT calculation, directly for the GW calculations 
as in Ref. ITU. Spin off-diagonal elements in the Green 
function and the self-energy are thus fully taken into ac- 
count. Finally, we investigate the behavior of the valence 
and conduction bands not only on certain lines of the BZ, 
but also in the whole symmetric mirror plane. 

The calculations and ARPES results are compared in 
Fig. [TJ The cut in Fig. [TJa) corresponding to /cf ^ = 
A , i.e., normal emission, has two features that dis- 
perse along kz. These are derived from the CB and VB. 
The CB is partially occupied by the degenerate bulk dop- 
ing of the sample. The data set also shows three features 
with non-dispersive two-dimensional (2D) character. The 
first one is the topological state at the Dirac point which 
is located at a binding energy of Eb ^ 0.45 eV. The sec- 
ond and third features are the so-called M-shaped state at 
Eb ^ 0.75 eV and the parabolic free-electron-like state 
at Eb ^ 0.25 eV (the latter is indicated in Fig. [TJa) 
as 2DEG - two-dimensional electron gas). These states 
have been interpreted as the 2D states formed by the 
quantization of the valence and conduction band states, 
respectively, in a potential well formed by the downward 
bending of the bands at the surface HHIii. Alter- 
natively, in Refs. 26l-[29| the emergence of the states is 
explained by an expansion of van der Waals spacings due 
to intercalation of surface-deposited atoms. The inter- 
calation of adsorbed Rb atoms, on the other hand, does 
not change the surface electronic structure significantly 
and so this point remains controversial [30| . The photoe- 
mission intensity of the surface-related features shows a 
resonant enhancement at the photon energies for which 
bulk states with a similar wave function periodicity along 
kz is observed. The topological state and the 2D state 
in the conduction band are enhanced at F whereas the 
M-shaped state is enhanced near Z. 

The best agreement between experimental and theoret- 
ical data is reached in the case of the LDA-fGVK calcu- 
lations with the relaxed atomic positions. However, the 
experimental band-gap value (0.332 ±0.055 eV) is not so 
well reproduced as in the LDA+GW case with experi- 
mental atomic positions (see Table H]). As is clearly seen 
in Fig. [B^a), the shown LDA results, which reflect the 
situation when the VBM is located in the mirror plane 
of the BZ (see Fig.[2j), are quite far from the experimen- 
tally observed valence-band edges. Also, it appears that 
all calculations overestimate the total width of the upper 
valence band. 

Similar to the conventional semiconductor systems, the 
GW corrections to the LDA band structure have mainly 
"moved" the conduction band away from the valence 



4 




r -0.1 -0.2 -0.3 -0.4 -0.5 r -0.1 -0.2 -0.3 -0.4 -0.5 

k- (A"') k- (A"') 

FIG. 3: (Color online) Projections of the lowest conduction 
band and the uppermost valence band in the mirror plane on 
the V — M direction of the two-dimensional BZ. The presented 
results are obtained for the experimental (upper row) and 
relaxed (lower row) atomic positions without (LDA) and with 
(LDA+GVK) the GW corrections to the LDA band structure. 



band on the energy scale. This fact is visually repre- 
sented in Fig. [21 where one can catch sight of only slightly 
changed contours of the bands in the mirror plane upon 
including the many-body corrections. Nevertheless, there 
is a crucial point that drastically distinguishes the TI 
from the conventional semiconductors. This point is the 
band inversion near the center of the BZ. Owing to this 
inversion, as was recently shown in Ref. 8, the mentioned 
movement apart decreases the "penetration" of the bands 
into each other near the F point, i.e., the hybridization 
due to spin-orbit coupling is reduced, and, as a result, 
both the band- inversion region in /c-space and the F-point 
band gap become smaller. 

For bismuth selenide, the mentioned mechanism of the 
many-body corrections finally leads to shifting the VBM 
to the F point, as is clearly seen in Fig.[3l where an effect 
caused by the relaxation of the atomic positions is also 
shown. On the DFT level, the relaxation results in a shift 
of the VBM along the mirror plane towards the center of 
the BZ and a reduction of the indirect band gap. With 
the inclusion of the GW corrections, the relaxation causes 
a smaller direct band-gap value and a larger upward shift 
of the conduction band beyond the band-inversion region. 

Fig.lHshows the GGA band structure without and with 
the GW corrections. As compared to the LDA-based 
calculations, the presented data allow one to address the 
question of how different approximations to the XC func- 
tional affect the GW results. It is worth noting that we 
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FIG. 4: (Color online) Same as in Fig.[3l but without (GGA) 
and with (QGA+GW) the GW corrections to the GGA band 
structure. 



TABLE I: The calculated band-gap values in eV for Bi2Se3. 



XC 


type 


unrelaxed 


relaxed 


LDA 


indirect 


0.30 


0.25 


LDA+Giy 


direct 


0.34 


0.19 


GGA 


indirect 


0.31 


0.28 


GGK+GW 


direct 


0.30 


0.21 



have obtained not only similar behaviour of the bands un- 
der study for both reference band structures (LDA and 
GGA), but even qualitatively close GW^-results on the 
valence-band edge, as one can judge from Fig.[TJ Besides, 
the performed DFT+GW^ calculations with different XC 
functionals and relaxed atomic positions give practically 
the same band-gap value (see Table H]). Also, the RPA 
dielectric constant £oo turns out to be quite stable upon 
changing the XC functional and relaxed atomic positions 
and is consistent with experimental data (see Table SI of 
Supplemental Material [23]). 

In conclusion, we have presented ARPES measure- 
ments for bismuth selenide which are aimed at a study 
of the behaviour of the bulk valence-band in the high- 
symmetry mirror plane in the bulk BZ. We have probed 
a large range of /c-space by taking data at a dense mesh 
of angles and photon energies (only data along f — M 
shown here). We have theoretically considered the effect 
of many-body corrections within the GW approximation 
on the band structure, in particular, on the lowest con- 
duction and highest valence band, and demonstrated how 
ah initio GW results depend on the DFT reference band 
structure and on atomic positions of Bi and Se atoms in 
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the rhomb ohedral crystal structure of Bi2Se3. As a net 
result, we have arrived at theoretical and experimental 
data which are in good agreement and indicate consis- 
tently that Bi2Se3 is a direct-gap semiconductor with the 
VBM located at the T point. 
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